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Abstract 


Energy  storage  devices  such  as  composite  rotors  are  built  wifli  radial  precompression  to 
enhance  mechanical  performance;  however,  the  preload  might  decrease  due  to  the  viscoelastic 
behavior  of  materials  at  elevated  temperatures.  In  this  investigation,  an  analytical  solution  is 
developed  to  study  the  viscoelastic  problem  of  thick-walled  cylinders.  The  analysis  accounts  for 
ply-by-ply  variations  of  properties,  fiber  orieutations,  and  temperature  gradients  througb  die 
thickness  of  cylinders.  Fiber-reinforced  composite  materials  generally  illustrate  extreme 
anisotropy  in  viscoelastic  behavior.  The  viscoelasticity  exists  mainly  in  matrix-dominant 
properties  such  as  transverse  and  shear  while  die  fiber-dominant  properties  behave  more  like 
elastic  mediums.  Accordingly,  the  viscoelastic  characteristics  of  composite  cylinders  are  quite 
different  from  those  of  isotropic  cylinders.  Currentiy,  finite-element  packages  such  as 
ABAQUS,  ANSYS,  and  DYNA3D  are  not  suitable  for  the  viscoelastic  analysis  of  composite 
cylinders  because  of  the  lack  of  anisotropic  viscoelastic  elemrats.  The  prestress  in  the 
hoop-wound  fibers,  which  generates  radial  compression  in  the  cylinders,  might  decrease  due  to 
Poisson's  effect  alone  from  the  creep  behavior  in  the  transverse  properties  of  composite.  The 
result  also  shows  the  effects  of  layup  constraction  and  fiber  orientations  on  the  anisotropic 
behavior  of  composite  cylinders. 
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1.  Introduction 


Conq)Osite  mataials  are  currently  used  for  highly  efficient  rotors,  as  an  energy  storage  device. 
For  these  applications,  prestresses  are  built  in  during  fabrication  of  the  cylinders  through  a  press-fit 
procedure  to  enhance  the  mechanical  performance.  In  a  high-speed  rotor,  centrifugal  force  resulting 
from  the  rotation  during  operation  generates  tensile  stresses  in  the  radial  and  circumferential 
directions.  Since  the  con^osite  rotors  are  primarily  reinforced  in  a  circumferential  direction,  the 
radial  tensile  is  critical  to  the  ultimate  performance  of  the  rotors  due  to  the  weakness  of  radial 
strength.  Accordingty,  it  is  crucial  to  build  and  maintain  a  certain  level  of  radial  con5)ression  in  the 
rotor  prior  to  operation.  Reinforced  con^site  materials  are  also  used  for  construction  of  lightweight 
gun  barrels.  Conposhe  gun  tubes  are  generally  built  with  a  metal  liner  and  a  con:q)osite  overwrap. 
Because  of  a  thick- walled  construction,  there  is  a  stress  gradient  through  the  thickness  of  barrels  as 
subjected  to  propellant  pressure.  For  a  gun  barrel  application,  there  is  a  stress  gradient  through  the 
wall  thickness.  The  stress  gradient  has  to  be  minimized  to  achieve  the  optimal  structural 
performance.  Con^osite  gun  barrels  must  be  prestressed;  therefore,  the  con5)osite  overwrap  is  in 
a  state  of  tension  to  induce  conq)ression  in  the  steel  liner. 

Polymer  materials  in  general  behave  viscoelastically,  which  means  that  creep  associated  with 
stress  relaxation  occurs  in  a  long  period  of  time,  especially  at  elevated  ten^eratures  (Ferry  1961). 
The  associated  stress  relaxation  in  the  composite  wall  results  in  loss  of  the  prestresses  and  potentially 
leads  to  structural  failures  or  downgrade  in  performance.  Rber-reinforced  con5)osites  with  polymer 
matrix  carry  the  same  characteristic,  but  behave  more  con^licatedly  due  to  anisotropy.  For 
carbon-fiber  conq}osite,  fiber-dominant  properties  such  as  modulus  in  the  fiber  direction  are  sort  of 
elastic.  The  viscoelasticity  exists  mainfy  in  matrix-dominant  properties  such  as  transverse  moduli  and 
shear  propaties.  Because  of  the  extreme  anisotropic  characteristic,  a  conqwsite  cylinder  will  behave 
quite  differently  fi'om  one  built  with  isotropic  materials.  The  objective  of  this  investigation  is  to 
develop  an  analytical  method  to  study  the  viscoelastic  behavior  of  the  thick-walled  conq)osite 
cylinders.  The  analysis  can  be  applied  to  the  design  of  flywheel  machinery  and  overwrapped 
con:q)osite  gun  tubes. 
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To  date,  activities  in  the  research  of  viscoelasticity  have  mainly  concerned  isotropic  materials 
included  in  the  studies  by  Muki  and  Sternberg  (1961),  Schapery  (1964),  Williams  (1964),  and 
Christensen  (1982).  These  basic  theories  of  viscoelasticity  were  then  extended  io  the  area  of 
heterogeneous  and  anisotropic  materials  for  a  variety  of  applications.  Hashin  (1965)  used  the 
effective  relaxation  moduli  and  creep  con^liances  to  define  the  macroscopic  viscoelastic  behavior  of 
linftar  viscoelastic  heterogoieous  media  and  its  inplementation  in  viscoelastic  modeling.  The  general 
formulation  of  linear  viscoelastic  boundary  value  problems  of  composite  materials,  including  the 
thermal  viscoelastic  problems  for  thermorheologically  single  materials  and  the  applications  of  the 
correspondence  princ^te,  were  examined  by  Schapery  (1967).  Rogers  and  Lee  (1964)  investigated 
the  viscoelastic  behavior  of  an  isotropic  cylinder. 

In  the  following  research,  the  linear  quasi-static  viscoelastic  behavior  of  a  thick  laminated 
composite  cylind^  with  an  elevated  tenperature  change  is  studied.  The  analysis  accounts  for 
pty-by-pfy  variation  of  properties,  tenq)erature  changes,  and  fiber  orientations.  The  thick  cylinder  is 
assumed  to  be  in  the  absence  of  thermomechanical  coupling  and  to  be  in  the  state  of  generalized  plane 
strain,  such  that  all  the  stress  and  strain  con^onents  are  independent  of  the  axial  coordinate  (Tzeng 
and  Chien  1994).  Moreover,  due  to  the  nature  of  axisymmetry,  all  the  stress  and  strain  conponents 
are  also  independent  of  the  circumferential  coordinate.  The  mechanical  responses  of  this  thick 
composite  cylinder  will,  thoiefore,  onfy  have  to  satisfy  the  governing  equation  in  the  radial  direction. 
Invoking  the  Boltzmann  superposition  integral  for  the  con:q)lete  spectrum  of  increments  of  anisotropic 
material  constants  with  respect  to  time,  the  thermoviscoelastic  constitutive  relations  of  the  anisotropic 
conposite  cylinder  can  be  derived  in  integral  forms.  Since  the  thick  conq)osite  cylinder  is  subjected 
to  a  constant  elevated  tenqjerature  change,  and  boundary  conditions  are  all  independent  of  time, 
formulations  of  the  linear  thermal  viscoelastic  problem  can  have  forms  identical  to  those  of  the 
corresponding  linear  thermoelastic  problem  by  taking  advantage  of  the  elastic-viscoelastic 
correspondence  principle.  In  other  words,  all  of  these  integral  constitutive  equations  reduce  to 
algebraic  relations  that  are  very  similar  to  those  developed  for  thermoelastic  media  when  they  are 
Laplace  transformed  by  means  of  the  rule  for  convolution  integrals.  The  thermoelastic  analysis  can 
thus  be  used  to  derive  the  transformed  thermal  viscoelastic  solutions  in  the  frequency  domain. 
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2.  Viscoelastic  Formulation 


Consid^  a  filament-wound  axisymmetric  thick  con^site  cylinder  consisting  of  N  layers  with  the 
axial  coordinate  z,  the  radial  coordinate  r,  and  the  circuniferential  coordinate  0,  as  shown  in  Figure  1. 
This  conposite  cylinder  has  the  inner  radius  a,  the  outer  radius  b,  and  the  length  L.  The  Boltzmann 
superposition  integral  of  a  stress  Oy  (i,  j  - 1, 2, 3)  and  strain  Cy  (i,  j  -  1, 2, 3)  relation  for  an  isothermal 
viscoelastic  problem  with  a  constant  ten:5)erature  increase  AT  and  the  thermal  expansion  coefficient 
a^is 

o,j(t)  =  -  P,j(T,t)AT  ,  (1) 


where  Cy*  (T,t)  is  the  relaxation  modulus  dependent  on  temperature  T  and  time  t. 


qf(T(t),t)  =  C^*(To,X(t)) 


X(t)  = 


t 

ar  (T(t))  ‘ 


Here,  To  is  the  base  ten^erature,  and  a^  is  the  ten^erature-shift  factor.  PyCT.t)  is  given  by 
By(T,t)  -  Qj“(T,t)  •  Ou .  It  is  often  desirable  to  use  the  inverse  form  of  the  constitutive  relation  (1), 

€,.(t)  =  AiJ'(T,t  -  +  ilrjj(T,t)AT  ,  (2) 


kl 

where  i|fij(T,t)  is  the  tensor  product  of  the  creep  compliance  Ajj  (T,t)  and  the  thermal  creep 
coefficiait  Since  the  elevated  tenq)erature  change  AT  is  constant  above  some  reference  value 
in  time,  the  relaxation  moduli  and  creq)  con^liances  are  evaluated  at  that  reference  ten^jerature 
regardless  of  whether  or  not  the  material  is  thermorheologicalfy  single  by  eirq)loying  the  temperature- 
shift  factor. 
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Figure  1.  Continuity  Cylindrical  Coordinate  System  (r,  q,  z)  of  a  Laminated  Cylinder 
Containing  N  Layers. 


The  Laplace  transfonn  of  a  function  f[t)  is  defined  as 


f  =  f(s)  =  r  e-“f(t)dt, 
Jo 


(3) 


where  s  is  the  Laplace  transform  variable.  Applying  (3)  with  the  convolution  rule  to  (1)  and  (2) 
reduces  the  integral  constitutive  equations  to  the  following  algebraic  relations: 


and  the  inverse  form 


respectively,  where 


y.  _  pkl  ~  ft 

~  Tkl  _  r  AT 

e..  =  A  -  a..  +  li;..— , 

ij  u  ij  5  ’ 


(4) 

(5) 


Cj'  =  sC^/,  aJ*  =  sAJ*,  a,.  =  sa..,ij.  =  s^..,p„  =  Cj*-  o..,  and  =  aJ*  • 
Furthermore,  it  can  be  shown  that  Ja^J  =  JCy  j  * . 
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PnnsiHftr  a  corresponding  thrainoelastic  problem  with  the  transformed  displacement  con5)onents 
u,  V,  and  ^  in  the  axial  direction,  the  circumferential  direction,  and  the  radial  direction, 
respectively,  in  each  layer.  The  axisymmetric  character  of  the  thick  composite  cylinder  along  with 
the  assun:q)tion  of  the  state  of  g^ieralized  plane  strain  leads  to  a  sin^lified  displacement  field,  which 
reflects  the  circumferential  independence  and  only-radial  dependence  of  w , 

u(r,0,2)  =  u(r,z) 
v(r,0,z)  =  V(r,z) 
w(r,0,z)  =  w(r). 

Substituting  these  transformed  displacement  con5)onents  into  the  strain-displacement  relations 
and  invoking  the  compatibility  conditions,  one  can  derive  explicit  forms  of  u  and  v .  Since  each  layer 
of  the  thick  laminated  cylindra*  is  cylindricalfy  monoclinic  in  reference  to  the  global  coordinates,  there 
is  no  coupling  between  transverse  shears  and  other  deformations.  Accordingly,  the  vanishing  ^ear 
traction  boundary  conditions  and  interface  continuity  conditions  generate  zero  out-of-plane  shear 
tractions  and  shear  strains  for  each  layer.  Moreover,  owing  to  the  absence  of  torsional  deformation, 
the  transformed  displacement  components  u  and  V  become 

u  =  S®z 

V  =  0,  (7) 

where  the  constant  quantity  E  ®  has  the  physical  interpretation  of  transformed  axial  strain  of  a  layer. 
In  feet,  S  ® ,  according  to  the  present  formulation,  also  represents  the  transformed  axial  strain  of  the 
entire  conqwsite  cylindw.  The  calculation  of  S  ®  requires  the  knowledge  of  end  boundary  conditions 
and  will  be  given  later.  Likewise,  solving  for  w(r)  requires  the  information  of  transformed  strain 
conqponents,  the  constitutive  equations,  as  weU  as  the  equilibrium  equations. 

The  previously  transformed  displacement  field  gives  the  transformed  strain  con^onents  in 
cylindrical  coordinates 
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^  dMr) 
"  dr 


^ee  " 


_  m 


“  dz 


Cft,  =  €„  =  6  =  0 

Or  zr  z 


The  unabridged  form  of  the  constitutive  equation  (4)  for  each  layer  in  cylindrical  coordinates  with 
the  radial  coordinate  r  normal  to  the  plane  of  symmetry  is  expressed  as 


^11  ,9i2  Ss 

§12  §22  §23 

^13  ^23  ^33 


0  0 
0  0 
0  0 


0  0  0  -44  ^45  0 


0  0  j3 

^16  ^26  ^36 


"55  0 


0  0 


Furthermore,  from  the  previous  discussions,  it  can  be  shown  that  two  of  the  three  equilibrium 
equations  are  satisfied  automatically.  The  only  nontrivial  equilibrium  equation  is  the  one  in  the  radial 
direction: 


-  Oee 


Substituting  (6),  (7),  and  (8)  into  (9),  the  transformed  stress  con:q)onents  and  dgg  are 
obtained  in  terms  of  the  transformed  radial  displacement  w.  Incorporating  the  resulting 
and  000  functions  with  (10)  gives  a  nonhomogeneous  Euler  differential  equation  of  wfor  a  layer 


^  ^  f  (B„-  B„)  -  (C„-  C.,)*-  , 

dr^  dr  C  s 
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where 


Solving  (11)  for  wyields 


= 


W 


=  Ajr^ 


A  *» 


V 


where 


Wp  =  fjE-r  +  f^r 


(12) 

(13) 


(14) 


and  Aj  and  A^are  coefficients  to  be  determined  from  boundary  and  continuity  conditions. 

Finally,  it  is  understood  that  the  initial  condition  of  the  original  thermoviscoelastic  problem  is 
displacement-free  state  of  rest.  The  boundary  conditions  are  of  free  tractions  and,  hence,  of  free 
transformed  tractions  on  both  inner  and  outer  circular  surfrces: 

®ir  =  ®er  =  ®  ^ 


On  both  end  surfaces,  stress  resultants  are  zero: 


N 


E  /'• 

k=i  '■ 


0„rdr  =  =  0  at  z  =  0,  L, 


(16) 
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where  Tj  and  r,,  are  inner  and  outer  radii,  respectively,  of  the  kth  layer.  The  continuity  conditions  at 
each  intflrface  between  two  adjacent  layers  require  continuous  radial  traction  and  continuous  radial 
displacement  at  any  instant  as  shown  in  Figure  2.  Thus,  when  written  in  the  transformed  form,  they 
become 


and 


ITyO 


(17) 

(18) 


where  k  =  1,  -,  N  -  1,  and  subscripts  i  and  o  denote  inner  and  outer  surfaces,  respectively. 

Accordingty,  the  formulation  accounts  for  ply-by-ply  variations  of  material  properties  and 
tenq)erature  change.  The  matrix  form  numerical  solution  procedure  with  parallel  con:q)utmg 
techniques  resolved  the  con:q>lexity  and  time-consuming  calculation  procedures  in  a  Laplace  transform 
of  a  multilayered  con^osite  cylinder  (Chien  and  Tzeng  1995). 


Ill  :  Radial  Displacement  o  :  outside  surface 

"on.:  Radial  Stress  i  :  inside  surface 


Figure  2.  Continuity  Boundary  Conditions  of  Radial  Displacement  and  Stresses  at  the 
Interfaces  of  Adjacent  Layers. 
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3.  Relaxation  of  Thermal  Stresses 


The  time-dependent  thermal  viscoelastic  behavior  of  a  100-layer  AS-4/3502  graphite  epoxy 
rnnnpnsite  cylinder  subjected  to  a  tenq)erature  increase  AT  -  150“  C  is  examined.  The  con^osite 
cylinder  has  an  innar  radius  a  ■■  3.5  in,  an  outer  radius  b  —  4.1  in,  and  a  thickness  of  each  layer 
h  ■■  6.0  •  10’^  in.  Stacking  sequence  is  given  as  [0/30/60/90]25  from  inside  out  with  the  0  direction 
coinciding  with  the  axis  of  the  cylinder.  The  creep  properties  of  AS-4/3502  graphite  epoxy  and 
AS-4/PEEK  thermoplastic  laminate  with  a  fiber  volume  fiaction  of  0.67  were  investigated  at  different 
tenperatures  by  Kim  and  Hartness  (1987).  Figure  3  shows  the  transverse  and  shear  con^liance  of 
the  composite  at  various  ten:5)eratures.  Increase  of  compliance  with  time  due  to  creep  behaviors  of 
material  was  found  at  elevated  temperatures. 


-0.5  0  05  1.0  1.5  2  0 


LOG  I,  MIN 


Figure  3.  Time-Dependent  Transverse  and  Shear  Compliances  at  Various  Temperatures  (Kim 
and  Hartness  1987). 


In  this  study,  the  least-squares  curve  fitting  was  used  to  express  the  creep  con^liances  fi’om  the 
original  AS-4/3502  data  at  T  -  150®  C  in  a  power  law  form 

$22(0  =  [1.7051  +  1]  S22  (19) 

and 

S^(t)  =  [11.3076  +  1]  S^.  (20) 

where 

S22  =  7.5328- lO'Vi 
=  1.3834- 10-®/psi- 

The  compliance  in  the  fiber  direction,  Sn  -  5.9  •  lO'Vpsi,  and  Poisson’s  ratios,  V12  -  V13  -  0.3, 
V23  -  0.36,  are  assumed  to  be  time  independent.  Thermal  expansion  coefQcients  of  the  con^osites 
in  three  princ^al  directions  are  an  -  -0.5  -  10'®/°C  and  a22  (tss  -  40.0  •  10‘®/®C,  where  the 
negative  value  indicates  shrinkage  with  tenqierature  increase. 

Figures  4  and  5  show  radial  displacement  and  radial  stress  profiles  across  the  thickness  of  the 
cylinder  at  three  instants,  t  -  10"^  min,  t  -  lO^'®  min,  and  t  -  10*^  min,  respectively.  The  radial 
tractions  and  displacemrats  satisfied  the  continuity  conditions  at  the  inter&ces  of  layers  at  all  instants. 
The  “saw”  shaped  radial  stress  distribution  results  firom  the  discontinuity  of  material  properties  and 
various  fiber  orientations.  The  radial  displacement,  w(t),  reaches  a  steady  state  at  t  -  10^^  min  due 
to  a  long-term  creep  behavior.  The  radial  displacement  of  most  layers  reaches  a  constant  value, 
except  at  the  innermost  and  outermost  portions  of  the  cylinder.  The  fi-ee  traction  boundary  at  the 
surface  of  cylinders  causes  the  gradients  in  the  radial  displacements.  A  similar  phenomena  is  also 
obsCTved  in  the  radial  stresses  showing  a  constant  level  of  stress  distribution  in  the  inner  portion  of 
cylinder  away  fiom  boundary  layers  at  t  -  10*^  min.  This  long-term  creep  characteristic  reflects  the 
proposed  power  law  form,  (19)  and  (20),  of  the  creep  con^liances. 
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Figure  4.  Radial  Displacement  Profiles  Across  the  Thickness  of  the  Cylinder  at  Three  Different 
Instants:  t  - 10*’  min,  t  -  lO*"*  min,  and  t  ■  lO*’  min. 


Figure  5.  Radial  Stress  Profiles  Across  the  Thickness  of  the  Cylinder  at  Three  Different 
Instants:  t  10*’  min,  t »  10^  min,  and  t  10^  min. 
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The  hoop  stress,  OgoO)  through  the  thickness  of  the  cylinder,  is  illustrated  at  three  instants  in 
Figures  6(a),  6(b),  and  6(c),  respectively.  There  exist  two  distinct  values  (discontinuity)  of  Oe9(t) 
across  each  interface  of  two  adjacent  layers  due  to  the  various  fiber  orientations  through  the 
thickness.  Thrae  is  a  trend  showing  the  creep  behavior  of  the  hoop  stress  profiles  fi’om  t  -  10'^  min 
to  t  -  10‘*  min.  The  hoop  stresses  in  60°  and  90°  layers  show  a  fairly  steep  gradient  across  the 
cylinder  thickness  initially.  However,  the  gradient  gradually  disappears  as  time  approaches  infinity. 
Likewise,  Figures  7(a),  7(b),  and  7(c)  present  the  axial  stress  o^(t)  profiles  in  the  radial  direction  at 
the  three  specified  instants,  respectively.  These  figures  exhibit  similar  creep  and  stress  relaxation 
behaviors  as  those  of  hoop  stresses  O09(t).  However,  the  axial  stresses  in  30°  and  60°  layers  show 
an  increase  outwardly  fi'om  the  inside  of  the  cylinder  in  the  early  stage.  Figiures  8(a),  8(b),  and  8(c) 
are  the  inplane  shear  stresses  Ojo(t)  through  the  thickness  of  the  cylinder  at  the  three  specified 
instants,  respective^.  The  iiq)lane  shear  stresses  of  0°  and  90°  layers  vanish  in  the  entire  time  history 
due  to  generalized  plane  strain  assun^tion.  A  drastic  change  occurs  in  the  30°  and  60°  layers, 
resulting  fi'om  the  combined  effects  of  Poisson’s  ratio  and  creep  characteristics  of  materials. 

4.  Relaxation  of  Mechanical  Stresses 

In  the  following  study,  a  calculation  was  performed  to  simulate  a  conqiosite  cylinder  pressed  fit 
on  a  rigid  steel  mandrel  as  shown  in  a  schematic  in  Figure  9.  The  0.6-in-thick  con^osite  cylinder  has 
a  layup  construction  of  [90/90/90/0]25, 75%  in  hoop  and  25%  in  axial  direction.  The  contact  pressiure 
is  1,000  psi  at  the  interfice  of  conposite  and  mandrel,  which  results  fiom  the  initial  interference.  The 
assembfy  is  thoi  heated  and  held  at  93°  C  constantly  for  creep  study.  In  modeling  of  the  experiment, 
the  inner  radius  of  the  conq)osite  overwrap  was  held  constantly  to  allow  creep  and  stress  relaxation 
of  the  cylinder  occurring  at  93°  C.  The  AS-4/epoxy  corrqiosite  material  properties  were  obtained 
fiom  Kim  and  Hartness  (1987).  A  power  equation  was  obtained  by  a  least-square  curve-fit  technique 
as  follows  (Tzeng  1997): 


S^(t)  =  7.553  ♦  10-''  (t)o.oi75  (21) 

and 

Sgg(t)  =  1.383  *  10-^  (t)®'*^ .  (22) 
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Radius  (inch) 


Figure  6.  Hoop  Stress  Profiles  Across  the  Thickness  of  the  Cylinder  at  Three  Different 
Instants:  (a)  t  - 10’^  min,  (b)  t  -  min,  and  (c)  t  «=  10“  min. 


Axial  stress  (psi)  Axiol  stress  (psl)  Axial  stress 


Figure  7.  Axial  Stress  Profiles  Across  the  Thickness  of  the  Cylinder  at  Three  Different 
Instants:  (a)  t  -  10'^  min,  (h)  t  -  min,  and  (c)  t  -  10^  min. 
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Figure  9.  Precompression  at  the  Inter&ce  of  a  Composite  Cylinder  and  Rigid  Mandrel  Can  Be 
Achieved  By  a  Press-Fit  Process  With  a  Proper  Interference. 

The  transv^^  and  shear  con^liances  are  in  a  unit  of  1/psi.  The  conq>liance  in  the  fiber  direction  is 
also  assumed  to  be  elastic,  as  discussed  in  the  previous  section. 

Figures  10  and  1 1  illustrate  the  creep  behavior  of  radial  displacement  and  relaxation  of  radial 
stresses  as  a  function  of  time,  respectively.  As  shown  in  the  plots,  a  significant  loss  of  preload  was 
in  the  conqwsite  cylinder  constructed  with  this  specific  material  and  layup,  subjected  to  an  elevated 
ten^rature.  Tests  have  been  poformed  using  similar  basis  graphite/epoxy  con^osite  materials,  and 
about  20%  of  preload  was  lost  in  a  few  weeks  at  elevated  ten:q)eratures.  In  general,  the  rate  of  loss 
of  preload  increases  as  the  ten^rature  increases.  The  stress  relaxation  in  hoop- wound  fiber  layers 
results  mainly  from  Poisson’s  effects  of  creq)  in  the  transverse  direction  since  the  properties  are 
assumed  to  be  elastic  in  the  fiber  direction.  The  stresses  will  not 
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RADIAL  DEFORMATION  @  93C 
(9Q/9Q/9Q/0) 
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Figure  10.  Creeping  With  Time  of  the  Radial  Displacement  Through  the  Thickness  of  a 
Composite  Overwrapped  Cylinder. 
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Figure  11.  Relaxation  With  Time  of  Radial  Stress  Through  the  Thickness  of  a  Composite 
Overwrapped  Cylinder. 
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Himinish  completely  due  to  elastic  properties  in  die  fiber  direction.  However,  the  reduction  and 
redistribution  of  hoop  stresses  will  cause  a  decrease  of  the  radial  compression  in  the  cylinders  in  a 
long  period.  As  the  radial  displacement  approaches  a  uniform  state  through  the  thickness,  the  radial 
stress  will  diminish  completely,  depending  upon  the  materid  properties  and  temperatures. 

The  viscoelastic  properties  are  commonly  not  available  for  most  composite  materials.  There  is 
a  need  to  develop  a  simple  and  accurate  estimation  of  the  viscoelastic  properties  of  composite 
material  for  the  analysis.  This  is  especially  important  fi»m  an  application  point  of  view  since  there 
are  numerous  composite  systems  available  through  a  combination  of  matrix  and  fiber.  The  following 
approach  is  proposed  to  provide  a  proper  method  of  material  selection  based  on  creep  and  stress 
relaxation  consideration.  The  creep  compliance  of  matrix  alone  can  be  measured  using  dynamic 
mechanical  analysis  (DMA).  Figure  12  shows  the  measured  data  of  ICl/Fibrite  977-2  resin  at  a 
reference  temperature  of  121  “  C.  The  data  can  tiien  be  curve-fit  with  a  time-dependent  function  as 
follows: 


SCO-S^xCt)®"^  (23) 

The  expression  basically  represents  the  time  dependency  of  the  matrix  material  of  composite 
material.  Since  die  viscoelastic  behavior  of  composite  material  is  mainly  dominated  by  the  matrix, 
it  is  reasonable  to  assume  the  transverse  and  shear  properties  of  graphite/epoxy  (977-2  resin) 
composite  following  the  power  law  as  equation  (23).  The  elastic  properties  that  are  generally 
available  can  then  be  substituted  into  the  expression  as  the  term  “S^”  to  obtain  the  viscoelastic 
properties  for  the  specific  composite  system.  The  proposed  approach  can  significantly  reduce 
engineering  efforts  in  material  testing  if  the  approximated  properties  can  be  applied  to  the 
viscoelastic  analysis  resulting  in  good  agreement  with  experiment 

In  the  following  study,  a  cylinder  analysis  was  performed  and  compared  to  experimental  data. 
A  cylinder  with  an  inner  diameter  of  10.0  in  and  a  thickness  of  1.2  in  was  fabricated  from 
graphite/epoxy  (Im7/977-2)  composite  material.  The  laminated  cylinder  has  a  ply  construction  of 
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Figure  12.  Creep  Compliance  of  Fibrite  977-2  Epoxy  Measured  by  DMA  at  a  Reference 
Temperature  of  121“  C  (Raw  Data  Are  Provided  by  CEM,  University  of  Texas  at 


Figure  13.  Relaxation  of  Interference  Pressure  of  a  10-in-Diameter  Composite  Cylinder  at 
Three  Time  Instants. 


an  interference  pressure  of  20  ksL  The  entire  assembly  is  then  kept  in  an  oven  at  121  ®  C.  The  hoop 
strain  at  the  inner  radius  of  the  mandrel  is  monitored,  which  provides  a  measure  of  the  preload 
variation.  A  loss  of  preload  was  found  to  be  18%  in  about  14  days.  The  result  of  analysis  is  shown 
in  figure  13,  which  indicates  a  fairly  good  agreement.  Although  the  validation  is  limited  to  a  single 
experiment  so  fer,  the  correlation  does  provide  confidence  of  inq)leinentation  of  the  proposed 
approach. 

5.  Conclusions 

An  analysis  has  been  developed  to  study  the  thermal  viscoelastic  behavior  of  thick-walled 
laminated  con^osite  cylinders  with  ply-by-ply  variation  of  anisotropic  viscoelastic  properties.  The 
relaxation  of  stresses  and  creq)  behaviors  is  properly  determined  in  a  cylinder  subjected  to  a  uniform 
temperature  change.  The  extreme  anisotropy  of  the  viscoelastic  behavior  of  composites  results  in 
very  investing  as  well  as  inportant  long-term  characteristics  of  cylinders,  which  cannot  be  modeled 
using  ciurentty  available  finite-elemoit  codes.  Creep  and  stress  relaxation  exists  in  the  fiber  direction 
even  though  the  fiber-dominant  properties  are  elastic.  This  is  mainly  due  to  the  contribution  of  the 
Poisson’s  effects  of  viscoelastic  transverse  and  shear  properties.  Furthermore,  the  analysis  is 
extended  to  study  the  viscoelastic  behavior  of  composite  overwrapped  cylinders  with  built-in 
mechanical  preloads.  The  relaxation  of  preload  in  the  composite  is  critical  to  the  service  life  cycle 
of  pplications  such  as  gun  tubes  and  composite  rotors,  which  are  designed  and  built  with  a  preload 
to  achieve  desired  mechanical  performance. 
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